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The equilibrium statistical mechanics of classical directed polymers in 2 dimensions is well known 
to be equivalent to the imaginary-time quantum dynamics of a 1+1-dimensional many-particle sys- 
tem, with polymer configurations corresponding to particle world-lines. This equivalence motivates 
the application of techniques originally designed for one-dimensional many-particle quantum systems 
to the exploration of many-polymer systems, as first recognized and exploited by P.-G. de Gennes 
[J. Chem. Phys. 48, 2257 (1968)]. In this low-dimensional setting interactions give rise to an emer- 
gent polymer fluid, and we examine how topological constraints on this polymer fluid (e.g., due to 
uncrossable pins or barriers) and their geometry give rise to strong, entropy-driven forces. In the 
limit of large polymer densities, in which a type of mean-field theory is accurate, we find that a 
point-like pin causes a divergent pile-up of polymer density on the high-density side of the pin and a 
zero-density region (or gap) of finite area on the low-density. In addition, we find that the force act- 
ing on a pin that is only mildly displaced from its equilibrium position is sub-Hookean, growing less 
than linearly with the displacement, and that the gap created by the pin also grows sublinearly with 
the displacement. By contrast, the forces acting between multiple pins separated along the direction 
preferred by the polymers are super-Hookean. These nonlinear responses result from effective long- 
ranged interactions between polymer segments, which emerge via short-ranged interactions between 
distant segments of long polymer strands. In the present paper, we focus on the case of an infinitely 
strong, repulsive contact interaction, which ensures that the polymers completely avoid one another. 
In a companion paper, we consider the effects of a wider set of inter-polymer interactions. 



I. INTRODUCTION 

The ensemble of configurations of a set of directed, one- 
dimensional objects embedded in (D + l)-dimensional 
space can be envisioned as the ensemble of worldlines of a 
corresponding set of quantum-mechanical point particles 
evolving in time in D spatial dimensions. In particu- 
lar, a standard mapping relates the classical canonical- 
ensemble equilibrium statistical mechanics of the set of 
directed one-dimensional objects to the imaginary-time 
evolution of the state of the corresponding set of point 
particles. This mapping was introduced and exploited 
by de Gennes^ in order to shed light on the equilibrium 
structure of directed fibrous polymers that are confined 
to two dimensions, thus providing a scheme for account- 
ing, nonperturbatively, for strong local polymer-polymer 
interactions that serve to prohibit configurations in which 
polymers cross. For de Gennes, this prohibition is ac- 
complished by asserting that the quantum particles are 
identical fermions (in his case free, noninteracting, and 
subject to periodic spatial boundary conditions), and are 
therefore subject to the Pauli exclusion principle. 

In a parallel development, a suite of powerful 
techniques — specifically, Bethe's Ansatz, bosonization, 
and quantum hydrodynamics — have been developed to 
address the quantum mechanics of one-dimensional sys- 
tems of many interacting particles or spins. The aim 
of the present paper and a companiorP one is to em- 
ploy these advances in quantum many-body (QMB) 
physics, together with the de Gennes analogy between 



the quantum-fluctuating, many-particle system and the 
classical, thermally fluctuating directed-polymer system, 
to uncover new information about the equilibrium struc- 
ture of systems comprising polymers that are either rigor- 
ously prohibited from passing through one another (i.e., 
noncrossing) or subject to other interactions, such as en- 
ergetic penalizations of crossings, or systems that allow 
for the presence of distinct species of polymers. In ad- 
dition, we apply these advances in technique to deter- 
mine the equilibrium forces acting on particles included 
in interacting polymer systems that serve to exclude the 
polymers from certain spatial regions, as well as the ef- 
fective forces that act between such particles as a result 
of their exclusion of polymers. A global theme of the 
present work is that, due to the reduced dimensional- 
ity of the polymer system, interactions dramatically in- 
fluence the structure and correlations that characterize 
the polymer system, and do so both, as we shall see, 
in topologically — and also geometrically — rich settings. 
This is a lesson already well known in the quantum- 
particle domain. Although we focus in this paper on 
polymer systems, our treatment applies to other two- 
dimensional statistical systems involving linelike degrees 
of freedom, such as wandering steps edges on crystal 
surfaces^, dynamically growing interfaces in the Kardar- 
Parisi-Zhang universality class^, and vortex lines in pla- 
nar type-II superconductors^. 

The present paper is organized as follows. In Sec.[ll]we 
introduce the two-dimensional interacting directed poly- 
mer system. We also describe its mapping to a one- 
dimensional quantum analog, and address the statistical 



2 



observables that may be derived via this mapping. In tion {a; n (-)}^L 1 of the polymer system: 



Sec. Ill we describe how to impose and analyze topologi- 
cal constraints on the polymer system, and we discuss in 
detail the resulting effects on the free energy and struc- 
ture of the polymer system. 



IV 



we summarize 



In Sec. 

our results and provide some conclusions. In a compan- 
ion papeil^, we consider more general classes of polymer 
interactions. For such systems, we analyze the resulting 
interpolymer correlations as well as the response of the 
polymer system to impurities such as free or fixed lines 
or particles, which are not topological in character. In 
addition, we discuss the application of the technique of 
bosonization, familiar from quantum many-body physics, 
as a tool for characterizing the universal behavior of di- 
rected polymer systems. 



II. DIRECTED 2D CLASSICAL EQUILIBRIUM 
POLYMERS AND EVOLVING ID QUANTUM 
PARTICLES 



Directed 2D classical polymers in thermal 
equilibrium 



We consider a system of N two-dimensional (2D) di- 
rected polymers, indexed by n = 1, . . . , N, that are non- 
crossing^. The configuration of the n th polymer is de- 
scribed by x n (t) , where r is the coordinate along the di- 
rected axis of the system (which we call the longitudinal 
direction) , and x n gives the location of the n th polymer 
in the perpendicular direction (which we call the lateral 
direction); see Fig. [I] We take the energy cost of the de- 
flections of the polymers from the longitudinal direction 
to be 



.4 



N 



2 ^ 



, jo 



n=l 



dr (d T x n ) 2 , 



(2.1) 



a N r L 

[{*«(•)}]= o £ / dr(d T x n (r)) 2 
n=l 

i N r L 

+ t£ / dr^(x n (T)) 
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+T £ / drV{x n {r)-x nl {T)), (2.2) 

in which we have also included an external (or one-body) 
potential term According to this model, in the ab- 
sence of polymer-polymer interactions or external poten- 
tials the polymer configurations have a thermal distri- 
bution that is Gaussian, in the sense that increments in 
their deflections [i.e., x n (r + St) — x n (r)] are indepen- 
dent Gaussian random variables having mean zero and 
variance SrkgT/A, where T is the system temperature 
and ks is Boltzmann's constant, which we generally set 
to unity via a suitable choice of units. We note that we 
shall not be considering the dynamics of the polymer sys- 
tem, so we do not need to take note of the kinetic energy 
of the polymer system. Additionally, we characterize the 
system via conditions on the configurations of the poly- 
mers at their ends, via the distributions P l ({x n }) and 
P' ({x n }), which respectively give the probability densi- 
ties for the configurations {x n } of the polymer ends at 
t = and t = L. 
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where L is the extent of system in the longitudinal di- 
rection, and A is the bending energy per unit length, 
which penalizes configurations for straying from the pre- 
ferred (i.e., longitudinal) direction. The polymers must 
be stiff in a way which we will define later to permit us 
to neglect higher-order terms in this expression. In ad- 
dition to the bending energy, we include an interaction 
V between the polymers, which we take to be trans- 
lationally and parity invariant, and sufficiently short- 
ranged that it may be taken to operate only between 
monomers (i.e., polymer segments) having common r co- 
ordinates. In fact, we shall often take the interaction to 
be purely local, in which case it would take the form 
V(x„(t) - x u '(t)) = c5{x n {r) - av(r)), where 5{x) is 
the one-dimensional Dirac delta function. Thus, we ar- 
rive at the following energy functional U of a configura- 



FIG. 1. The paths {x n (r)} describe a possible configuration 
of the directed polymer system. Thermal fluctuations permit 
the system to adopt energetically disfavored configurations. 
When polymers appear to intersect in the (x, r) plane, in 
reality one crosses over the other by exploiting the presence 
of a third dimension. 

We take the polymer system to be at thermal equilib- 
rium at inverse temperature (3. Thus, we have for the 
canonical ensemble partition function^ 

z[pf,p i } = J dixiypfaxi^dix^p'axi}) 

{ Xn {L)=Xl} 

x Jv[{x n (-)}]^ U[Mm , ( 2 -3) 
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which depends functionally on P l and P* . Here, 
the measures are defined via d{Xl}d{X^} = 

Uti dX[ dX], and V [{*„(•)}] = J[n=x V [f »(•)]• To 
complete the definition of this multiple path integral we 
also need to impose some form of lateral boundary con- 
ditions on the polymer configurations. We return to this 
point in Sec. |IID| 

As for the thermal expectation value (0[{x n (-)}}) of 
a generic observable (i.e., a functional of the polymer 
configuration) O [{x n (-)}], this is given by 

{o) = z[pf,PT 1 J d{xf}pf({xf})d{xi}p i ({xi}) 

{x n (L)=xl} 

© [{*«(•)}] e-^ x ^O[{x n (r)}}. (2.4) 

{x„(0)=X;} 

As we have noted above and shall see below, owing to the 
low dimensionality of this thermally fluctuating polymer 
system, even short-ranged interactions produce qualita- 
tive alterations of the structure and correlations that it 
exhibits, relative to those exhibited by its noninteracting 
counterpart. Moreover, even interactions that are weak, 
microscopically, are fundamentally nonperturbative, in 
that they induce correlations that are long-ranged. 



B. Mapping to ID quantum particles 

Let us turn now to the consideration of a one- 
dimensional quantum system of N nonrelativistic par- 
ticles each of mass m and having coordinates {g r n }^_ 1 , 
subject to a one-body interaction <£>(</„(£)) and to a 
(translationally and parity invariant) two-body inter- 
action V(q n — q n ')- For this system, and introduc- 
ing the (unsymmctrized) simultaneous particle-position 
eigenkets |{^}), a matrix element of the imaginary-time 
propagator ({Xl}\e~ Hr / h \{X l n }) can be expressed as the 
following Feynman integral over paths {q n (t)}; (see, e.g., 
RefP): 



({^}|e- 



HT/h \{K}) 

= J f [{««(•)}] e^^-)}]/", (2.5) 



where the Euclidean action Se is given by 
r j- n n 



s E = f dt{]T!(^„) 2 + £<%„) 

*'° n=l n=l 

+ Y, V(q n (t)-q n ,(t))}. (2.6) 



Kn<n'<JV 



The "non-Lagrangian" sign of the interaction term in 
Eq. ( |2.6[ ) and the terminology of the Euclidean action 
reflect the fact that we are considering imaginary-time 



propagation, in which the paths of the Feynman inte- 
gral are referred to as imaginary-time world-lines. As 
usual, the propagator can be used to construct the transi- 
tion amplitude between generic initial and final quantum 
states and 

(*/|e-^|«*) = J d{*>}d{**}<*/|{*>}) 

^{{X^e-^/^iKmX^m. (2.7) 

In order to relate this quantum-mechanical system to 
the classical polymer systems, we make the following 
identifications. We match the wave functions with the 
probability distributions, i.e., we choose 



P l ({Xn}), 
P f ({Xn})- 



(2.8a) 
(2.8b) 



Note that the identification is not between classical and 
quantal probability distributions but, rather, between 
classical probability distributions and quantal wave func- 
tions. Thus, to be appropriate, the wave functions should 
be restricted to being real, non-negative, and integrat- 
ing to unity. The equivalence between the quantal and 
classical problems is completed by adding the following 
identifications: 



T = h/3, 
t = Tt/L, 
q n (t) = x„(r), 
m = h 2 p 2 A/L. 



(2.9a) 
(2.9b) 
(2.9c) 
(2.9d) 



In order to maintain the analogy to the polymer system, 
we refer to the ID quantum system as having a width 
rather than a length, and denote the widths of both sys- 
tems by w. Then we have the result that is central to 
this paper, viz., that the quantal matrix element of the 
imaginary-time evolution operator is equal to the classi- 
cal partition function of the polymer system: 



(^f f \e- HT ^ i ) = Z[P f ,P l ] 



(2.10) 



In particular, the imaginary-time worldlines of the quan- 
tum particles correspond to the configurations of the di- 
rected polymers, and the quantum fluctuations of the 
particle system (the strength of which is governed by h) 
correspond to the thermal fluctuations of the polymer 
system (the strength of which is governed by 1//3). 

In order to apply many of the techniques of QMB 
physics it is useful to impose a choice of quantum statis- 
tics upon the initial and final quantum states and 
The condition that the initial and final wave 
functions ({A'„}| , I' t ) and ({J^^ } | \E f ^) be non-negative, 
so that they may match the polymer end distributions, 
Eqs. (2.8), precludes the direct choice of fermionic statis- 



tics (although, as we shall see, such statistics can be ap- 
plied indirectly, following a Jordan- Wigner-type trans- 
formation). Thus, we choose to consider situations in 
which the polymer endpoint distributions P % and P' 
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are each symmetric functions — respectively under the ex- 
change of initial endpoints amongst themselves and fi- 
nal endpoints amongst themselves — from which it follows 
that the initial and final quantum end states be symmet- 
ric under particle exchange, and hence that they describe 
identical bosons®. The quantum Hamiltonian that cor- 
responds to the Euclidean action ( |2.6[ ) is invariant un- 
der particle exchange (i.e., PHP^ = H for all pairwise 
particle exchange operators P), and therefore under the 
imaginary-time evolution described by Q~ Hr / h bosonic 
many-body states remain bosonic. 

One valuable notion made accessible and sharpened 
via the tools of one-dimensional quantum many-body 
systems (such as bosonization, quantum hydrodynamics, 
and Bethe-Ansatz-rooted methods) is that of the emer- 
gent directed-polymer liquid (cf. Kafri et al. 5 ). This state 
is a classical analog of the Luttinger-Tomonaga liquid, 
which can be exhibited by interacting QMB systems in 
one dimension, and is qualitatively distinct from, e.g., 
the Landau liquid state of many-fermion systems, i.e., a 
state that can be exhibited by such systems in higher 
dimensions. Thus, we shall see, e.g., that the man- 
ner in which density correlations decay spatially in the 
emergent directed-polymer liquid resembles the space- 
time decay of superfluid fluctuation correlations in one- 
dimensional interacting QMB systems. 

C. Eigenfunction expansion of the imaginary-time 
propagator 

Due to the association of polymer configurations with 
particle paths, discussed in the previous subsection, it 
is the path- integral (i.e., covariant) formulation of the 
quantization of the system that is the one most clearly as- 
sociated with the physical degrees of freedom of the fluc- 
tuating polymer system. However, the quantum mechan- 
ics of the same particle system can also be formulated in 
terms of the time-dependent Schrodinger equation (i.e., 
via canonical quantization). Thus, it is straightforward 
to ascertain that the propagator can be expressed in 
terms of the expansion 

({xi}\e- m / n \{xi}) 

over the exact normalized eigenfunctions {ipk} and cor- 
responding energy eigenvalues {Ek\, with many-body 
quantum numbers fc, of the associated many-body quan- 
tum Hamiltonian 

N 2 N 
n— 1 n— 1 

+ V{x n -x n ,). (2.13) 

l<n<n'<N 

An important special case is made evident via this ex- 
pansion. Suppose one is concerned with a statistical- 



mechanical expectation value involving polymer observ- 
ables all taken at a single value of the longitudinal co- 
ordinate r and, moreover, obeying (t,T — t) 3> h/AE, 
where AE is the spacing between the ground many-body 
state \tpgs) of the QMB system and its first excited state. 
In this case, it can be adequate to retain only the ground 
state in the eigenfunction expansion, i.e., to take 

«|^.)e- B "*/»tya (2.14) 

a situation referred to as ground-state dominance. In 
particular, within the ground-state dominance approxi- 
mation and far from the system ends, the equilibrium 
expectation value of the polymer density 

<£%-x„(t))) (2.15) 

n 

is given, as a function of lateral position x, by the 
quantum-mechanical ground-state expectation value of 
the density operator J2 n H x ~ Qn)> a result that holds re- 
gardless of the longitudinal boundary conditions on the 
polymer configurations. Thus, properties deep in the lon- 
gitudinal interior of a long system (i.e., one for which 
T ^> h/AE) are associated with the ground-state prop- 
erties of the quantum system. 

D. Noncrossing polymers 
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FIG. 2. Snapshot of a configuration for the case of noncross- 
ing polymers. Such polymers are not permitted to cross one 
another, but are otherwise noninteracting. 

Our focus is on systems comprising strictly noncross- 
ing polymers. Such systems do not adopt configurations 
in which any polymer crosses any other but, beyond this 
important element, they are not subject to any polymer- 
polymer interactions; see Fig. [2j As we shall see in a 
companion paper^l, QMB physics techniques enable the 
study of polymer systems having a wide range of inter- 
actions; we shall show there that in the presence of such 
interactions many of the results obtained in the present 
paper will continue to hold, at least qualitatively. To 
enforce noncrossing, the polymers are taken to feature 
an infinitely strong excluded- volume effect, so that the 
partition function contains only those paths for which 
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x n (r) 7^ s„'(r) for all n ^ n' . This restriction can be 
enforced for the polymer system via the inclusion of the 
interaction term 

V(x n (r) - av(r)) = c5(x n (r) - x n ,{r)), (2.16) 

with c/LA — > oo. The corresponding quantum system 
can therefore be taken to comprise many identical par- 
ticles, bosonic in their quantum statistics and subject 
to an interparticle interaction having the same form, 
Eq. ( |2.16[ ). This system is known as the hard-core (or 
impenetrable) point-like boson model. (A finitely strong 
interpolymer repulsion corresponds to a quantum Lieb- 
Liniger system^.) We mention that, despite its short 
range, this interaction is strong, and we may therefore 
expect the qualitative behavior of noncrossing polymer 
systems to be replicated in systems having more general 
interactions. 



1. From bosons to fermions 

It was demonstrated by Girardeau^ that any system of 
interacting, one-dimensional, bosonic particles for which 
the wave functions vanish when pairs of particles coin- 
cide in space can be mapped to an equivalent system 
of fermionic particles, subject to the same interaction 
V(xi — Xj) whenever Xj, ^ Xj. Girardeau's mapping is 
particularly useful in the case of hard-core boson sys- 
tems, as these can be mapped to systems of free fermions. 
Thus, as first noted and exploited by de Gennes^, the 
noncrossing condition on polymers can be accounted for 
entirely by the quantum statistics of fermions, without 
the need to include any interaction term. We mention 
that Girardeau's mapping preserves the modulus of the 
quantum wave function in the position basis; however, 
the forms of the bosonic and fermionic momentum-space 
wave functions are not preserved. This means that the 
local density of polymers is correctly described under the 
mapping but, e.g., quantities involving the slopes of poly- 
mer configurations — the analogs of the momenta of the 
quantum particles — are not. De Gennes applied this free 
fermion picture of two-dimensional noncrossing fibrous 
polymers to describe the structure of such systems. In 
particular, he showed that there is a logarithmic diver- 
gence in the "x-ray form factor" (i.e., the longitudinally- 
averaged correlator between lateral Fourier components 
of the density fluctuations of the polymer system), at a 
lengthscale associated with the Fermi momentum of the 
quantum system (or, equivalently, the mean lateral inter- 
polymer spacing). 



2. Ground state and ground-state dominance 



As noted in Sec. |II C the ground state of the quantum 
Hamiltonian plays a key role in the behavior of the poly- 
mer system over long distances. In order to make use of 
this idea, we now obtain the ground-state wave function 



for the many-hard-core boson system, subject to hard- 
wall boundary conditions, in a form that is convenient 
for the subsequent analysis. To do this, we begin with 
the ground-state wave function tpg S ({x n }) of a system of 
N hard-core bosons on a ring of circumference w, subject 
to periodic (rather than hard-wall) boundary conditions 
(see e.gP2); this is given by 



0iV(JV-l)/2 

n 



Kn<n'<N 



. 7T 

sin — ( x 

w 



(x„-av)|. (2.17) 



Such boundary conditions are appropriate for a system 
of polymers that lie on a cylindrical surface and are di- 
rected along the cylinder axis. Our aim, however, is to 
consider a system of polymers that are confined to a strip 
with hard-wall boundary conditions, and this system cor- 
responds to a quantum system also subject to hard-wall 
boundary conditions. To that end, we give the ground 
state wave function of a system of ./V bosons subject to 
vanishing boundary condtions at x n = ±w/2 (see Ap- 
pendix Ib| : 



1p gs ({x n }) = 



2 JV 2 /2 

i N / 2 y/W\ 



JV 

n 



cos 



n 



sin ■ 



sin ■ 



Kn<n' <N 



(2.18) 



The wave function ip gs reflects the inter-polymer repul- 
sion. Although the corresponding polymers are forbidden 
energetically only from actually intersecting one another, 
continuity and thermal fluctuations have the combined 
effect of "carving out" a spatial region around the poly- 
mers so that the probability of finding one polymer very 
near another (compared with the mean inter-polymer 
spacing w/N) vanishes as the square of the separation. 
A similar effect occurs near the hard boundaries, i.e., at 
x = ±w/2. The preceding results pertain to infinitely 
strong contact interactions. However, as is known from 
the work of Lieb and Linigep^, the physical properties of 
a system of bosons subject even to weak contact interac- 
tions differ nonperturbativcly from those of a system of 
free bosons. 

Whereas the lateral correlations amongst the poly- 
mer segments depend additionally on the quantum- 
mechanical energy eigenfunctions, the thermodynamic 
properties of the polymer system are determined solely 
by the spectrum of energy eigenvalues. The ground-state 
energy and energy spacing to the first excited state of the 
quantum system are, respectively, given by 



^9$ o 9 
o mw 



N = 



L 



AE 



2 V KT 

i -N 



6 w 2 P 2 A 
o L 



iV' 5 



w 2 /3 2 A 



N 



(2.19a) 
(2.19b) 



G 



for N ^> 1. The partition function of a long system is 
dominated by the term exp {—E gs T/h). 

Thus, the free energy density of a long system of non- 
crossing polymers is given, to leading order, by 



wL 



7T~ 



N 



1 



/3 2 A' 



(2.20) 



It is straightforward to show that, in contrast, the free 
energy density of a system of noninteracting polymers 
would be smaller by a factor of 3/N 2 . The noncrossing 
conditions essentially restricts each polymer to a region 
of width of order w/N, leading to a strong reduction in 
entropy and therefore increase in free energy. Even a 
weak inter-polymer repulsion would suffice to generate a 
free energy proportional to N 3 rather than N. 

The ground state energy is proportionate to the mean 
square polymer slope ((d T x n ) 2 ). This term must be small 
in order to justify our expression for the deflection energy 
which necessitates 



A N 

^ > — ■ 
1 w 



(2.21) 



We conclude this section on ground state dominance by 
remarking that the lengthscale over which the ground- 
state dominance approximation holds is t > L//3AE, 
or 



t > 



w 2 /3A 
N 



(2.22) 



III. TOPOLOGICALLY CONSTRAINED 
SYSTEMS OF POLYMERS 

A. Effects of a single pin on noncrossing polymers 

1. Constraints on the partition function due to a single pin 
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FIG. 3. A topological obstruction (a pin), located at (a;, r) = 
(x p ,Tp). Thermal fluctuations cannot carry a polymer across 
this pin, so a fixed number of polymers Nl pass to one side 
of the pin, and the remainder Nr = N — Nl pass to the other 
side. 

So far, we have been considering systems of polymers 
that are, from the geometrical and topological stand- 
points, trivial. We now turn our attention to systems 



that are subject to a pin; as we shall see, this renders 
them nontrivial, topologically. By a pin we mean a region 
of the polymer system, sharply localized near the point 
(x, t) = (x p ,Tp), at which the one-body potential expe- 
rienced by any polymer segment is taken to be so large 
and repulsive that the polymers never cross it during the 
course of an experiment. Thus, the pin serves as a topo- 
logical constraint and, because the polymers are directed, 
it partitions the configuration space of the polymer sys- 
tem into sectors labeled by the number N L of polymers 
that have the property that as they pass through the 
line t = t p they obey —w/2 < x(t p ) < x p ; see Fig. [3] 
Then, the corresponding number of polymers that pass 
the pin on its other side [i.e., obey x p < x(t p ) < w/2] 
is given by Nr = N — N L . We note that on the line 
t — T p the mean polymer densities to the left and right 
of the pin are, respectively, p L — Nl/{{w/2) + x p ) and 
pn = Nr/((w/2) — x p ). Evidently, the constraint cre- 
ated by the pin eliminates polymer configurations from 
the thermal ensemble, and thereby reduces the entropy 
of the system. For a generically located pin, the only 
configurations remaining in the ensemble correspond to 
large fluctuations of the original system. The pin there- 
fore raises the free energy of the system and, as a result, 
there is generically an equilibrium force on the pin. 

To determine the increase in the free energy due to 
the presence of the pin, it is convenient to analyze the 
partition function of the polymer system, restricted to 
having N L polymers constrained to pass on one side of 
the pin (as described more precisely in the previous para- 
graph), normalized by the unrestricted partition func- 
tion. This amounts to computing Eq. \2A\ with the ob- 
servable 0[{x n (-)}] given by 5{n l , J2n=i 6{xp-x n {Tp)\ , 
where 6(N, N') is the Kronecker delta function (i.e., 1 for 
N = N' and for N ^ TV') and &(•) is the usual Heaviside 
step function, which takes the value and 1, respectively, 
for negative and positive arguments. We now set about 
computing this free energy increase, as well as the impact 
of the pin on the spatial variation of the polymer density. 



2. How the constraint is reflected in the quantum 
propagator 



To compute various consequences of the presence of 
the pin we employ the mapping from the polymers to 
quantum-mechanical particles. This necessitates that we 
focus on the propagator, interrupted at imaginary time 
tp at which time we require that there be precisely N L 
particles in the range —w/2 < x < x p but with their 
exact locations left otherwise unspecified, i.e., 



d{x n }(^ f \c 



-H(T-t p )/h 



\{Xn}) 



x({x n }\e- m ^ i ). (3.1) 
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Here, the symbol C indicates the following constraints on 
the range of integration: 

—w/2 < xi, . . . ,xn l < x p 

< x Nl+ i, . . . ,x N < to/2. (3.2) 

In quantum-mechanical language, we are therefore to 
compute two factors: the amplitude to have precisely 
Nl particles at specific positions to the left of the pin 
(and the complement of particles at specific positions to 
the right of the pin) at imaginary time t p , given that the 
system was in the state l^) at imaginary time 0; and 
the amplitude to have the state |^/) at imaginary time 
T —t p , given that precisely Nl particles were at specific 
positions to the left of the pin (and the complement of 
particles at specific positions to the right of the pin) at 
imaginary time 0. 



Next, we invoke ground-state dominance, discussed in 



Sec. IIC which is justified by the restriction to values of 



the time t p that are far from the initial and final times 



(i.e., and T). Thus, the amplitude of Eq. (3.1), now 



normalized so as to yield the quotient of partition func- 
tions discussed in Sec. ITT1 becomes 



d{x n } (V ga \{x n }) ({x n }\V gs ) 



= J c d{x n }\({x n }\y gs )\ ; (3.3) 

the exponential factor exp {—E gs T /ft) has canceled be- 
tween the numerator and denominator, and \^ gs ) is the 
normalized ground state of the quantum system, the nor- 
malization being given by {9 gs \^/g S ) = 1. 



3. Transformation to the density profile 



To perform the constrained integration we invoke the form of the ground-state wave function, Eq. (2.18), and thus 
rewrite Eq. (3.3) as (choosing now units so that w = 7r) 



d{x n } \({x n }\^ g 



ir N Nl 



N 



d{x n } exp I ^^lncos 2 x„+ In ( sin x n — sin x T: 



l<n<n'<N 



(3.4) 



This (generically high-dimensional) integral is conveniently analyzed by exchanging the discrete particle coordinates 
{x n } for the continuum density function p, via 



1 N 

P({ x n}) = — ^2S(X-X n ), 



(3.5) 



and replacing the multiple integration over polymer coordinates by functional integration over the configurations of 
the polymer density, via which we obtain 



2 N 

d{x n }\{{x n }\^ gs )\ 2 = ^ w ^^ 



r/2 



N 2 



tt/2 



Vp exp I N I dx p(x)hicos x+—— dx dx'p(x) p{x') In 



sin x — sin x 



(3.6) 



-7T/2 -7r/2 

Here, the symbol C indicates the condition on p that corresponds to the pin constraint C, viz. 

dxp(x) = N L /N, 

tt/2 

but also the normalization and positivity conditions, which are respectively given by 

r n/2 



tt/2 



dx p(x) = 1; 

p(x) > 0, for - tt/2 < x < tt/2; 



(3.7a) 

(3.7b) 
(3.7c) 



both of which follow from the definition of p as a density, Eq. (3.5). Thus, in Eq. (3.6) we have arrived at a 
convenient formulation of the partition function in the presence of the pin constraint. The convenience results from 
the replacement of particle coordinates by the collective density field, which enables us to employ the techniques 
of the calculus of variations, and thus to determine the value of the density that yields the dominant (in the sense 
of Laplace's method) contribution to the functional integral. A more rigorous transition from individual polymer 
coordinates to a density function via the method of collective coordinates-^ would yield an additional 0(N In N) term 



in the exponential in Eq. (3.6). Our current level of accuracy suffices when all length scales, such as the displacement 



of the pin from its equilibrium position, are long compared to the average inter-polymer spacing. 

It should be noted that in Eq. (3.6) we have an effective long-ranged (logarithmic) repulsion between polymer 
segments, which replaces the original, local interaction. Two polymers that are laterally separated at t = t p will 
nevertheless interact indirectly with one another, owing to the nonzero probability that they would encounter one 
another at some distant value of r. Thus, as we are effectively integrating out all polymer degrees of freedom that do 
not lie on the line r = t„, we are generating an additional inter-polymer repulsion. We have thus mapped our two- 
dimensional system with its short-ranged interactions on to a one-dimensional system having long-ranged interactions. 
Such a description is adequate for addressing large fluctuation phenomena on the line t = t p that arise from a pin or 
other constraints. 

I 



4- Obtaining the density profile 



The factors of N in the exponent of Eq. (3.6) sug- 



gest that we may treat this functional integral via the 
functional version of Laplace's method, and therefore via 
finding the maximum of the leading term of the expo- 
nent. To this end, we approximate the partition function 



max exp 
c 



(3.8) 



where the free energy functional describing the effect of 
the pin AJ- [p] is defined via collecting the leading terms 
in Eq. 



AF[p] = -N 2 T Tn2 

7r/2 

N 2 T f r _ 2 

dxdxp(x)p(x) In (since— sin x) 



(3.9) 



-7I-/2 

We note that this expression for the free energy scales 
with N as iV 2 . As they do for the overall free energy of 
the system, the inter-polymer interactions qualitatively 
change the response of the polymers to the pin. 

Our next task is determine the approximate value of 
the mean density profile p{x) that minimizes the free en- 



ergy, subject to the constraints discu ssed in Sec . HI A 3 
Upon implementing the constraints ( 3.7a|3.7b ) via La- 



grange's method of undetermined multipliers, we arrive 
at the following first-order stationarity condition: 



/ dx p(x) In I sin a; — sini| 

J-w/2 

+Ai + A 2 6(x p - x) = 0. (3.10) 



It is important to note that this first-order condition is 
only required to hold at values of x for which p{x) is 
greater than zero, owing to the positivity condition (3.7c) 



that p is required to obey. In the absence of the pin con- 
straint, the A2 term would be absent, and the condition 
would be trivially satisfied via a uniform polymer density 
profile. This uniform density profile is also the correct 
mean density profile in the presence of a pin for the ex- 
ceptional case in which the pin demands that the average 
polymer densities on the two sides of it be equal to one 



another (i.e., for the case in which pl = PR = N/w. We 
refer to such a constraint as an equilibrium pin. Generi- 
cally, however, the pin has a profound effect on the mean 
density profile. 



To analyze the stationarity condition, Eq. (3.10), it is 



convenient to make the following transformations of the 
dependent and independent variables: 

p(x) — > Q(s) = p(x)/ cosx, (3.11a) 
x — > s = sin ir, (3.11b) 

so that dx = ds/y/l — s 2 , and similarly for x. In terms of 
these new variables the stationarity condition becomes 

J ds Q(s) In \ s — s\ 

+Ai + A 2 6(s p - s) = 0, (3.12) 

where s p = sin x p . This form of the stationarity condition 
must be met at values of s for which Q(s) > 0. 

If we were to ignore the positivity condition that a 
physically acceptable density profile Q must satisfy then 
it would be straightforward to solve the integral equa- 
tion (3.12) for all s (by following a technique that we 
shall, in fact, eventually adopt; see, e.g., RefP^). How- 
ever, the result we would obtain for Q would be invalid, as 
it would diverge to negative infinity near the pin. Thus, 
we search for a solution that violates the first-order con- 
ditions over one or more segments of — 1 < s < 1; within 
these segments Q(s) =0 and we refer to any such a seg- 
ment as a gap. 

Nonequilibrium pins necessarily cause compression of 
the the polymers either on one side of the pin or on the 
other. Without loss of generality, we may take the com- 
pressed region to correspond to — 1 < s < s p ; then the 
rarified region corresponds to s p < s < 1. For the pro- 



file that minimizes the free-energy functional (3.9), it is 



physically reasonable to assert that any gap would form 
on the rarified side and, furthermore, that it would lie di- 
rectly adjacent to the pin and extend over a region that 
more strongly disfavors the presence of polymers, as a re- 
sult of its proximity to the dense region of polymers that 
is created by the pin on the other side of the pin. (These 
assertions will be verified as we proceed.) Thus, we look 
for non-negative profiles for which we require Q(s) = 
for s lying within the range s p < s < s g . If a gap were to 
be present anywhere else, the system free energy could be 
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lowered by moving polymer density into that gap. Thus, 
any solution found that (i) has this gap and (ii) elsewhere 
satisfies the constraints on p(x) is the unique minimize! - 
of our effective free energy. 
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FIG. 4. Transverse variation of the equilibrium directed 
polymer density along the line r = r p for a system subject to a 
pinning constraint located at (x, r) = (—to/4, t p ), the location 
of which is marked by a black dot. The physical polymer 
density is the plotted quantity scaled by Ntt/w. Within the 
present approximation scheme, a finite gap (i.e., an area of 
zero polymer density) extends from the pin. The dashed line 
represents the equilibrium polymer density in the absence of 
a pin. 

As we shall see, the formation of a single gap yields a 
suitable profile (i.e., one that satisfies the positivity con- 
straint automatically). What remains, then, for the case 
of a single pin, is to determine the value of the density 
profile and, in particular, the width of the gap (s g — s p ) 
that together ensure that the necessary conditions on the 
profile are obeyed by it. Following the approach reviewed 
in RefP^, we have that the family of solutions to this in- 
tegral equation is given in terms of parameters Aq and 

M by 



(A Q + Ais)sgn(s - s g ) 
(> ( s) = I ^(1- s 2 ){s- s p )(s- s g ) ' 

0, 



for — I < s < s r , 



or s g < s < 1; 



for s p < s < s g . 



Transforming back to polymer coordinates, using 
Eq. ( |3.11 1, we thus have p(x) given by 



(Aq + A\ sina;) sgn(sinx — sina; 9 ) for — tt/2 < x < x p 
*/ (sinx — sin x p ) (sin x — sinx g ) ' or x g < x < tt/2; 



0, 



for x p < x < x g . 



The next step in determining p(x) is to adjust A and 
Ai to ensure that there is no divergence at x = x g ; this 
requires that Ao + Ai sinx g = 0. Physically, this choice is 
motivated by the expectation that the equilibrium den- 
sity does not diverge on the rarified side of the pin. Next, 
we invoke the normalization condition (3.7b) and thus 



determine that A\ — 1/tt; and, finally, we adjust x g to 
ensure that the pin constraint (3.7a I is met. Thus, we 



arrive at the profile that dominates the constrained par- 
tition function: 



p{x) 




for — 7r/2 < x < x p 
or Xg < x < tt/2; 



for x p < x < Xg. 



Restoring the physical lengths, this becomes p(x) 



1 / sin(7ra;/w) — sin(7nr fl /u?) for — w/2 < x < x p 
w y sin(7ra;/w) — sm(TTx p /w) ' or x g < x < w/2; 



.0, 



for x p < x < x g . 



Note the essential features of this solution for the polymer 
density along the line r = t p , as shown in Fig. [4] through- 
out the width of the system (i.e., for —w/2 < x < w/2) 
the density is non-negative; at the system edges (i.e., 
x = ±w/2) the density is finite; at the pin (i.e., for 
x = x p ) the density has a square- root divergence ap- 
proaching x p from the compressed side and a square- 
root vanishing approaching x g from the rarified side; and 
within the segment x p < x < x g the density is zero. It 
is striking that merely as a result of inter-polymer in- 
teractions that are local in the two-dimensional plane 
the topological restriction presented by the pin causes 
the opening up of a finite gap in the polymer density 
and, in particular, that the reach of its impact extends 
over many times the intrinsic inter-polymer separation, 
at least at the level of the present mean-field type of 
approximation. The mechanism responsible for this is 
that — over a longitudinal distance that is nonzero — the 
polymers are energetically disfavored from entirely fill- 
ing in the gap by the cost in bending energy they would 
have to incur to depart from their long-distance equilib- 
rium positions. Note that while the mean density p(x), 
as calculated within our approach, is zero within the gap, 
fluctuations can only increase it, leading to there being 
a small positive polymer density within this region. 

A closed-form expression for the gap edge x g associated 
with a partitioning may be found in terms of elliptic in- 

For a pin near its equilibrium 
the pin displacement x p — x f = 6 and the 



tegrals, as in Eq. (Bla 
position x 



gap x p — Xg = 7 are both small. Then the gap size is 
related to pin displacement to leading order by 



In 



2w 



cos 



w 



+ 0( 7 ). (3.13) 



Thus for a pin near its equilibrium position the gap size 
is sublinear in the pin displacement. 

In Fig. [5] we show the gap size more generally as a 
function of the pin position and the parameter v = 
(pL — Pr)/(pl + Pr), which we introduce to characterize 
the imbalance between the densities on either side of the 
pin. When all the polymers are on one side of the pin 
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FIG. 5. Dependence of the width of the gap that opens in 
the polymer density adjacent to a pin on the pin location and 
polymer density imbalance. Note that the gap width increases 
from zero as the pin position is varied, parametrically, from its 
equilibrium position. At the equilibrium position (p^ = pr), 
the gap vanishes and the polymer density is spatially uniform, 
as it is for a system in the absence of a pin. 



(in which case v = ±1) then the gap size is simply the 
pin displacement. More generally, the polymers on the 
rarefied side expand towards the pin so that the gap is 
smaller than the pin displacement. 



B. Force on a pin 

1. Calculation of the force 

Having determined the mean density profile on the 
lateral line through the pin, viz. p(x, r)| T=T , we now 
return to our formulation of the dominant contribu- 
tion to the partition function, and hence the free en- 
ergy, to determine the increase in the free energy due 
to the pin, viz. AJ 7 . To do this, we begin with 



the free energy functional Eq. (3.6) and seek to com- 
pute its value at the mean density profile. To sim- 
plify this computation, we employ the first-order sta- 



tionarity condition (3.10) to eliminate the combination 

— J*_^y/ 2 dx' p(x') In \$in(irx/w) — sm(nx' /w)\ in favor of 

Ai + A2 0(x p — x). We also use the stationarity condition 
evaluated at x = ±w/2 to obtain the Lagrange multipli- 
ers Ai and A2 in terms of the (known) mean profile. Then, 



in the resulting expression for AJ 7 , we use the normal- 
ization of p, Eq. (3.7b), and the constraint on it that the 
pin introduces, Eq. (|3.7a|), and thus arrive at the result 



AT(x p ,N L ,N,T,w) 



-N 2 T 



w/2 



w/2 



dx p(^x . x p ^ x g J 



:((N L /N) In 2(1 + sin(7ra/u;)) 

+(N R /N) In 2(1 - sm(nx/w))), (3.14) 



We remind the reader that the gap edge location x g is not 
an independent variable but is determined in terms of the 
independent variables, via the pin constraint, Eq. (3.7a). 



As a special case, we first consider the situation in 
which all polymers lie to one side of the pin, i.e. Nl = N, 
so that the pin can in effect be taken to be a septum 
emerging normally from one wall of the system. In this 
case, the formula for the free energy increase, Eq. (3.14), 
simplifies to the explicit form 



AT - 



N 2 T 



dx 



s'm(Ttx/w) — 1 



-w/2 



sin(TTx/w) — sin(TTX p /w) 
x In 2(l + sin(-Kx/w)) 



= -N 2 T\n ( [1 + sm(irx p /w)] /2) . (3.15) 

Remaining with the case Nl = N, we note that for mild 
compressions (i.e, those obeying (w/2) — x p <C w) the 
force exerted on the pin is Hookean in nature: 



AT- 



-NT 



{w/2) - x p 



(3.16) 



as established by making Taylor expansion of Eq. (3.15) 
about x = w/2. On the other hand, for strong compres- 
sions [for which (w/2) + x p -C w], so that the polymers 
are forced through an opening the width of which is only 
a small fraction of the full width of the system, we have 
the form 



AF = 2N 2 T In ( - . 

7T V (w/2) + X 



(3.17) 



Residing, as it does, beyond the linear-response regime, it 
is not surprising that this form is non-Hookean. Indeed, 
a similar term, with N^ replacing A^ 2 , dominates the free 
energy for highly compressed systems with polymers on 
both sides of the pin. 

In the more general case, in which Nr > polymers 
pass to the right of the pin, the force on the pin obeys 
(see appendix |B| 



dT 



ttN 2 T ( sin(TTX p /w) — sin (irx g /w) 



iu 



cos (nXp/w) 



.(3.18) 



For a pin near its equilibrium position this force is simply 



proportionate to the gap size x p 



and the above 
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FIG. 6. Increase in free energy as a function of pin position, 
for various values of the number of directed polymers Nl 
passing to one side of the pin. The physical free energy is the 
plotted quantity scaled by N 2 T. The free energy cost due to 
the pin is less than harmonic for a generic pin that is near its 
equilibrium position, and it diverges logarithmically as the 
pin nears an edge of the system. The curves plotted here 
are for pins having equilibrium (and free energy-minimizing) 
positions — Aw, — .3w, — .2w ... Aw. 



expression may be integrated to obtain the leading term 
in the free energy 



4 V w ) 



2w 



(3.19) 



Note that although the free energy of the pin near its 
equilibrium position grows faster than quadratically in 
the gap size 7 it is sub-Hookean in (i.e., grows slower 
than quadratically with)the pin displacement x p — x e . 
In contrast, in the limit in which the pin position ap- 
proaches the boundary of the system and highly com- 
presses the polymers, the free energy diverges logarith- 
mically. Fig. [6] shows the free energy as a function of 
pin position for polymer partitionings corresponding to 
various equilibrium positions of the pin. 



C. Effects of a barrier on noncrossing polymers 

The approach we have developed so far addresses the 
case of a topological constraint created by an infinitcs- 
imally wide pin. We now generalize the approach to 
allow consideration of the case in which the constraint 
is a barrier of finite width, so that polymers are pro- 
hibited from passing through an extended line segment 
(x,t)\ t =t p (with— w/2 < < x < 2:^ < w/2), as shown 
in Fig. [7J We note that the effect of single barrier is, from 
the standpoint of directed-polymer statistical mechanics, 
entirely equivalent to the effect of a pair of pins, pro- 
vided that polymers are forbidden from passing between 
the pins. 

To analyze the situation of a barrier and, in particu- 
lar, to compute the change in free energy arising from 




(-w/2,0) 



(0,0) 



(w/2,0) 



FIG. 7. A barrier (i.e., a laterally oriented topological ob- 
struction of nonzero width) that constrains the number of 
directed polymers that pass on either side. By using the tech- 
niques developed in the present paper one can readily analyze 
the effects of such barriers. 



the presence of this barrier, we return to the task of min- 
imizing the free energy (3.9), but we replace the pin con- 
straints ( 3.7a|3.7b ), by the following ones, appropriate 
for a barrier in which precisely Nl polymers pass to the 
left of the barrier and Nr pass to its right: 



— 7r/2 
7T/2 



dxp{x) = N L /N, 



dxp(x) = N R /N. 



(3.20a) 
(3.20b) 



These constraints along with normalization and positiv- 

x L 

ity imply the barrier condition: J x £ dx p(x) = 0. 

Physically, it is evident that there are two distinct sit- 
uations. Consider a barrier of given width, and with a 
partitioning specified by Nl and Nr. Relative to the sit- 
uation without the barrier, the polymers are compressed 
on at least one side of the barrier, and possibly both. 
Let us focus on a compressed side, at which the polymer 
density diverges, and imagine shrinking the barrier into 
a pin located at Xp. On the other side of the pin, there 
would now be a gap in the polymer density, the width 
of which is determined, as before, by the pin constraints. 
Now, imagine widening the pin into a barrier. As long 
as the barrier width does not exceed the gap width, the 
polymer density profile would not change, remaining at 
the density profile associated with a single pin at Xp. In 
effect, the barrier resides in the gap created by the pin, 
so the fact that the barrier width is finite has no impact. 

The second situation follows when the barrier width 
does exceed the width of the gap created by the pin. Now, 
both ends of the barrier are in contact with polymers. 
Repeating the integral-equation analysis of Sec. |III A4| 
we find the resulting density profile (choosing units as 
before so that w = n) to be given by p(x) = 

(An + Ai sinx) sgn(sinx - sinXp ) for - tt/2 < x < x^ 
(sinx — sin 2;^) (sin 2; - s'mxp 1 ) or x p < x < ir/2: 



0. 



for x^ < x < Xp . 
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FIG. 8. Equilibrium directed-polymer density along the 
line containing a barrier. The physical polymer density is 
the plotted quantity scaled by Ntt/w. The spatial extent 
of the barrier is indicated by the thick line. Note that for 
the chosen value of the polymer partitioning, the polymer 
density diverges on both sides of the barrier. The dashed line 
represents the density of the polymer system in the absence 
of the barrier. 



The two constants are now set by the numbers of poly- 
mers passing on either side of the barrier. The area of the 
"gap" is now precisely the area excluded by the barrier, 
and the polymer density diverges on either side of it, as 
shown in Fig. [Hj 

Although, in general, the effect of the barrier on the 
free energy is more complicated than the effect of a single 
pin, there is one case that may be addressed analytically. 
Consider an extended barrier of width i< to, the mid- 
point of which would be an equilibrium pin if b were set 
to 0. In this case, the equilibrium polymer density on 
either side of the barrier resembles the polymer density 
around a single pin that has been displaced a distance b/2 
from an equilibrium position x — w/2. The free energy 
cost of such a small barrier is then given by 



AT 



7T 2 (N 

8 I w 



(3.21) 



D. 



Effects of multiple pins and/or barriers on 
noncrossing polymers 



We now have all the tools we need to address multiple 
pins and even multiple finitely-wide barriers, as long as 
we continue to restrict these obstacles to lying on a com- 
mon line r = T p . Consider, then, M pins at the ordered 
locations —{w/2) < x\ < X2 < ■ ■ ■ < xm < (w/2), such 
that Ni polymers are constrained to pass between pin I 
(located at xg) and the nearest pin (or wall) to its left 
(i.e., at smaller x). N\ and Nm+i respectively denote 
the number of polymers passing to the left of the left- 
most pin and the right of the rightmost pin. (We remind 
the reader that any barrier of finite width may be treated 



as an adjacent pair of pins with no polymers passing be- 
tween them.) For such situations, the cost in free energy 
is determined, as usual, by maximizing the logarithm of 
the appropriate partition function Z, given by (choosing 
here units so that w = ir) 



In Z ~ J d> x J dx p(x) p(x) In ( sin x 



sin x 



(3.22) 



over the density profile p(-), subject to the constraints 
imposed by the pins and/or barriers. The symbol C now 
indicates the following collection of constraints for all I: 



c> 



dxp(x) = Ng/N, 
p(x) > 0, for 



(3.23) 

(tt/2) < x < (tt/2), (3.24) 



where Cg indicates that the integration range runs to the 
I th obstacle from the obstacle or wall that precedes it (or, 
in the rightmost case, from the rightmost obstacle to the 
right wall). The barrier constraints are implemented via 
a collection of Lagrange multiplier terms, which augment 
In Z and are given by 



M+l 

E 



Xp 



dx p(x) 



Nt 
N 



(3.25) 



As with the cases already treated in Sec. |III A 4[ func- 
tional differentiation — of the free energy functional aug- 
mented by the Lagrange multiplier terms — yields an in- 
tegral equation that is solvable for a non-negative, and 
therefore physically acceptable, p(-), provided we (i) al- 
low for the possibility of gaps in the polymer density 
profile, and (ii) determine their necessity and location 
by implementing the constraints on the numbers of poly- 
mers passing between the various obstacles. In general, 
some pins will lie within the interior of gaps; other pins 
will have a gap form on one side. 

For cases involving more than one or two pins, the 
process of determining the gap structure that permits all 
constraints to be satisfied becomes quite tedious, but it 
should always yield a unique result for the density pro- 
file that minimizes the free energy. Moreover, denoting 
by {pj} the sets of points at which the polymer density 
diverges and by {gk} the set where it increases gradually 
from zero density, we may continue as before and thus 
obtain the mean polymer density in regions outside the 
gaps as 



p(x) 



_1 Y[ k [sm(Trx/w) ~ sin(7rg fc /w)] 
w y Y[ Asin(nx/w) — sin(npj/w)] 



(3.26) 



E. Connection to quantum hydrodynamics 

In describing the effect of a single pin, we have focused 
specifically on the equilibrium structure of the polymer 
system on the line r = r p that passes though the pin. To 
describe the polymer behavior away from this line (i.e., 
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for t 7^ t p ) is equivalent to describing the imaginary- 
time evolution of the quantum system away from this line 
that is consistent with the density profile p(x) at t = r p . 
As this profile differs substantially from the equilibrium 
density profile for the pin-less case (i.e., a uniform density 
profile) the equivalent quantum system can be regarded 
as having undergone a large quantum fluctuation. 

Systems of nonlinear equations analogous to hydrody- 
namical equations have been used to describe the evolu- 
tion of one-dimensional systems of interacting particles 
around large fluctuations^. In terms of the particle den- 
sity and velocity fields, p(x, t) and v(x, t), these equations 
are 

dtp + d x (pv) = 0, (3.27a) 
d t v + vd x v = m~ 1 d x d p (p£(p)), (3.27b) 

where £(p) is the ground-state energy per particle, ex- 
pressed as a function of the density p. 
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FIG. 9. Variation of the equilibrium directed-polymer den- 
sity over the (x,t) plane in the vicinity of a pin constraint. 
The physical polymer density is the plotted quantity scaled by 
Nn/w. Note that the gap (i.e., the region in which the present 
approximation scheme gives zero for the polymer density) is 
roughly triangular in shape, as shown by the unshaded area 
of the inset. The divergence in the polymer density occurs 
only at precisely the location of the pin, depicted by a black 
dot in the inset. 

For a system of free fermions £(p) = (ti 2 it 2 / 6m) p 2 , and 
the hydrodynamical equations may be expressed as the 
complex Hopf equation 



dtw — iwd x w 
w(x, t) 



0, 



(3.28a) 

p(x,t) +iv(x,t). (3.28b) 



This equation can be solved numerically to describe the 
polymer density in the x-r plane resulting from a pin, 
with the mapping between quantum and polymer pa- 
rameters as in Eq. (2.9). The previously obtained mean 



density profile serves as a boundary condition for the 
field w. w(x,t p ) = (Hir/m) p{x) [v(x,t p ) — due to 



symmetry considerations] . As the hydrodynamical equa- 
tions are obtained by ignoring contributions associated 
with certain higher-gradient terms, they are expected to 
be quantitatively incorrect in the vicinity of any pins, 
because there the polymer density diverges (according 
to the mean-field type of approach that constitutes the 
bulk of the present paper) , and it is higher-order terms 
that dominate and control such divergences. Indeed, for 
a broad class of systems, nonlinearities can lead to shock- 
wave behavior even for initially fairly smooth density pro- 
files^, and even more readily for the divergent behavior 
around a pin. The polymer density resulting from a pin 
at x = — 0.06w with the partitioning Nl — .74N is shown 
in Fig. [9] We note that the divergence in the polymer 
density at the pin immediately subsides but the gap per- 
sists over a finite range in the longitudinal direction. At 
large longitudinal distances, however, the polymer den- 
sity profile returns to the uniform distribution of the pin- 
less system. 

A qualitative analysis of large quantum fluctuations 
indicates that the probability P(R) of a large fluctuation 
over a lengthscale R has the form P(R) ~ exp(— ai? 2 )P^. 
This corresponds to our finding in the present paper that 
when the pin displacement is small (compared to the 
overall system size) and all polymers lie on one side of 
the pin the free-energy cost is quadratic in the displace- 
ment. However, when the finite temperature of such a 
quantum system is taken into account, the probability of 
a fluctuation becomes P(R) ~ exp(— jR). The analog 
of a finite temperature in the quantum system is a finite 
length in the polymer system. For a system in which the 
pin coordinate r p is located near enough to another sys- 
tem feature, such as another pin or an end of the system, 
the ground-state dominance approximation fails, and one 
instead finds that the free-energy cost of the pin would in- 
crease linearly with the displacement of the pin from its 
equilibrium position. Thus, a longitudinally short sys- 
tem, or one having longitudinally distributed pins, can 
display a super-Hookean response when a pin is displaced 
a small amount from its equilibrium position. 

We remark that although the polymer system we have 
been exploring is formally equivalent to its quantum ana- 
log, the polymer system is more readily controllable. 
Large quantum fluctuations are, due to their rarity, diffi- 
cult to observe. In contrast, the probability of occurrence 
of the equivalent large thermal fluctuations of the poly- 
mer system can be measured indirectly, via the entropic 
force on a pin, which has the useful effect of forcing the 
system to assume what would otherwise be rare configu- 
rations. 



IV. CONCLUDING REMARKS 

Our central aim with this paper has been to develop a 
statistical-mechanical treatment of macroscopic systems 
of directed, two-dimensional, classical polymers in ther- 
mal equilibrium in the limit of infinitely strong repulsion 
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between the polymers. More specifically, our focus has 
been on the consequences of topological constraints — 
either point-like pins or spatially extended lines — that 
place limitations on the microscopic configurations ac- 
cessible to the polymer systems. These constraints have 
the effect of altering the equilibrium value of the spa- 
tial profile of the local polymer density and reducing the 
entropy (and, perforce, increase the free energy) of di- 
rected polymer systems and, when generically located, 
they give rise to forces that act on the constraints (and 
between them, if there are several of them) . 

The main technique that we employ to compute these 
effects is a mapping between two-dimensional systems 
of classical directed polymers in thermal equilibrium and 
the imaginary-time evolution of one-dimensional systems 
of quantal point-like particles in the ground state. (The 
particles can be taken to be either non-interacting iden- 
tical fermions or identical bosons subject to infinitely 
repulsive short-ranged interactions.) The determination 
of the alteration of statistical weight resulting from the 
presence of constraints for the polymer system is then 
ascertained via consideration of a rather unusual quan- 
tal amplitude: the matrix element of the imaginary-time 
evolution operator between the many-particle ground 
state and an appropriately partitioned configuration of 
the particle system. This technique readily yields^ re- 
sults for situations in which the constraints are located 
collinearly on a line that runs perpendicular to the direc- 
tion preferred by the polymers. By supplementing this 
technique with a quantum hydrodynamic approach we 
are able to establish the form of the polymer density not 
only along the aforementioned line on which the con- 
straints lie but also at points away from this line. 

Strikingly, in the limit of large polymer densities, in 
which fluctuations are small and a type of mean-field the- 
ory is accurate, we find that the effect of a point-like pin 
is to cause a divergent pile-up of the polymer density 
on the high-density side of the pin and a zero-density 
region (or gap) on the low-density side of the pin, the 
latter giving way to a continuously rising density beyond 
the gap. In addition, via the quantum hydrodynamic 
approach we find that the gap opened by the pin has 
a nonzero extent in the preferred direction, only gradu- 
ally narrowing and closing. (Fluctuation corrections to 
this mean-field theory are expected to replace the zero- 
density region by a region of very small polymer-density.) 
We also find for pins that are only mildly displaced from 
a zero-force (i.e., equilibrium) position the force acting 
on them is sub-Hookean, growing less than linearly with 
the displacement via a factor logarithmic in the displace- 
ment, and the gap created by them similarly grows sub- 
linearly with the displacement. By contrast, for multiple 
pins that are separated from one another along the direc- 
tion preferred by the polymers we find free energy costs 
that are super-Hookean, growing more than quadratically 
with the displacement. One can regard these nonlinear 
responses as resulting from the effectively long-ranged 
interactions between polymer segments that emerge via 



short-ranged interactions between long polymer strands 
in regions that reach far from the segments in question. 

In a forthcoming companion paper^l we shall discuss 
developments arising from the present work that encom- 
pass, inter alia, the impact of external potentials of vari- 
ous forms, the consequences of inter-polymer interactions 
that are finitely (as opposed to infinitely) strong or ex- 
tended in range, the behavior of distinct polymer species 
within host polymer systems, and and the application of 
powerful techniques from quantum many-body physics, 
including bosonization and the Bethe Ansatz. 
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Appendix A: N fermion ground-state wave function 

Consider a system of N noninteracting fermions mov- 
ing freely one dimension on the line segment < x < 7r 
and subject to homogeneous (i.e., homogeneous Dirich- 
let) boundary conditions at x = and x = ir. The nor- 
malized wave functions associated with the single-particle 
energy eigenstates of such a system are given by 



(j}j{x) = V2/7T sin(jx), j = 1,2,. 



(Al) 



In terms of these, the normalized A^-particle ground-state 
wave function of the A-fermion system may be expressed 
in terms of a Slater determinant, as follows: 



ip(xi, . . .,x N ) 



1 



det [<f>j(x k )], 



(A2) 



where [4>j{xk)] is the NxN matrix having (jk) element 
(f>j(xk)- Omitting N factors of y/2/n, the matrix [(j)j(xk)] 
takes the form 



/ sin^i 
sin 2xi 



sin X2 
sin 2x2 



sin xn \ 
sin 2xn 



(A3) 



ys'mNxi s'mNx2 ■•• sin Nxn J 



To simplify the evaluation of the determinant of this ma- 
trix, we add linear combinations of the rows of the ma- 
trix to other rows, a procedure that leaves its value un- 
changed. Specifically, we make the replacements 

4>j(x k ) -t cf>j{x k ) + <f>j_ 2 (Xk) - 2 cos cc at (f>j-l(x k ), 
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where we make the definition 4>j{xk) = for j < 1. The 
resulting matrix is then given by 



sinxi • • • sinxjv-i 


sinxjy 


2(cosxfc — cosxn) 







x sin((j - l)x fe ) 




o 



(A4) 



/ 



Note, specifically, the vanishing of the elements in the 
lower N — 1 rows in the iV th column. In evaluating the 
determinant of this matrix, we may extract N — 1 fac- 
tors of the form 2(cosxfe — cosx^v) from the rows in the 
lower left block of the matrix. Thus, we obtain the re- 
cursive result, linking Slater determinants for ./V-particle 
and (N — l)-particle systems: 



tetJM**)] =(-!) 



N-l 



JV-l 



(2cosxfe — 2cosxjv) x det [^-(ifc)]. (A5) 



k=l 



Next, we apply this relation recursively, and this arrive 
at the following, desired form for the TV-particle ground- 
state wave function: 



4>(x 1 ,X2 1 ■ ■ -,x N ) 



2N 2 /2 
TT N / 2 yfN~\ 

n ( c ° ss 

l<j<k<N 



N 
3 

COSXk 



If 



SIM, 



(A6) 



Note that for a system of hard-core bosons the corre- 
sponding ground-state wave function is simply the abso- 
lute value of this fermionic wave function. 



Appendix B: Elliptic integral representations of the 
pin constraint and free energy 

Given the form of the polymer density found in 
Sec. Ill A 4 the pin constraint Eq. ( 3.7a[ ) may be con- 
structed in terms of the elliptic integral of the third kind 
as 



Nr = 2(s g - s p ) 

N ny/(l- Sp )(l + S g ) 



xn 



1 + Sp (1 + S p )(l - Sg) 

1 + V (1 -*„)(! + a„) 



il(n, m) 



o l-ns 2 y/(l -ms 2 ) (I - s 2 ) 



(Bla) 



(Bib) 



We may describe how s g varies with s p for a given par- 
titioning by requiring that the partitioning given by Nr 
is constant with respect to variations of s p , so that 



ds„ 



8N R dN R ds g 



ds. 



ds g dsp 



= 0. 



(B2) 



By using Eq. (Bla), we thus find the result 



8Sg 



l + s g 



E 



(l+s p )(l-s g ) 
(l- Sp )(l+s g ) ; 

K ( a+s p )(i- Sg y 



(B3) 



where K and E are elliptic integrals of the first and sec- 
ond kind, respectively: 



K(k) EE 

£(fc)EE 



ds 



^{i-k 2 s 2 )(i-s 2 y 

1 , y/l - k 2 S 2 

as — . 

o VT^ 



(B4a) 
(B4b) 



For a pin infinitesimally displaced from its equilibrium 
position, so that s p — s g sa jcosx p with j <C 1, the 



right-hand side version of the pin constraint (3.7a I, may 
similarly be asymptotically expanded in 7 to obtain the 
condition 



N 



ds 



7 COS Xp 1 



111 



2 s — Sj 
2 cos 2 x p 



(B5) 



This equation yields the asymptotic relation between gap 
size and displacement given in Eq. (3.13). 



We turn our attention now to the free energy cost of 
the pin, Eq. (3.14), which can be written in terms of 
the scaled effective potentials </>£ and 4> R experienced by 
polymers on the left and right side of the pin, respec- 
tively, as 



N R 



N 2 T N^ N ™ 

cp L = J ds In [2(1 + s)]Q{s), 

4>r = J ds In [2(1 - s)} Q(s), 



(B6a) 
(B6b) 
(B6c) 



with the transformed coordinate s and polymer density 



Q(s) defined in the main text; see Eq. (3.11 1 . As we shall 
see, it is useful to represent (f>R as an elliptic integral, and 
to do this it is necessary to obtain a form without the 
logarithm. To this end, we extend s into the complex 
plane, and make use of the residue theory result 



<j> ds\n{\ - s)Q(s) = 0, 



(B7) 



where the keyhole contour T comprises the segments 
, T4 as shown in Fig. 10 Now, the integrand has 



the form 



16 




FIG. 10. The contour T used in Eq. (BIO I. Branch cuts run 
from s — — 1 to s — min(s p , s g ) and from s = max(s p , s g ) to 
s = +oo. The arc I\ has a radius e which will be taken to 
zero, and the arc F2 has a radius R which will be taken to 
infinity. 



ln(l -s)Q(s) = -(ln|l-.s|+iarg(l-.s) 

7T V 



(l_ s 2)( s _ Sp ) 



i<f>(s) 



4>(s) = ( arg(s - s g ) - arg(l - s) - arg(l + s) - arg(s - s p ))/2, 



(B8a) 
(B8b) 



which is analytic in the complex plane except on the branch cuts. Note that although the physical polymer density 
is zero in the gap, we are allowing Q(s) to be nonzero but imaginary in the gap, and complex in the complex s plane. 
We choose the branch cuts to run from s = —1 to s = min(s p , s g ) and from s — max(s p , s g ) to s = +oo, as shown in 
Fig.[T0} 

Whereas the integral along the inner contour T4 vanishes as its radius e goes to zero, the integral along the outer 
contour T2 (i.e., along |s| = R), does not vanish. In the limit R 3> 1 this contour integral becomes 



ds ln(l - s)Q(s) 



2tv 



(iRe i6 d9) kiR + i(6-w) 



1 



-2]nR. 



(B9) 



Thus, the contour integral of Eq. (|B10 1 yields 

2 1m U ~- / dsln|l - s\ Q(s + ie) 



ds 



ds 



In II - 



Q(s- 



ln|l - s 
-1 



Q(s + ie) 
ds In 1 1 — s\ Q(s — it). 



(BIO) 



Elementary cancellations among parts of these integrals (which is exact in this limit) yields 



occur, so that upon taking the limit R — > 00, Eq. (BIO 



p 00 


N 


/, * 






/"OG 











(s*-l)(s-s p ) 



(s 2 - l)(s + s p ) 



(Blla) 
(Bllb) 
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Using these representations of 4>l and 4>r, the free en- 
ergy may be expressed in terms of elliptic integrals, al- 
though the complete expression is rather complicated. 
Strikingly, however, a simple result follows for the force 
on t he pi n, —dTjdx v . In particular, by differentiating 
Eq. |B6a|) we find 



1 dT N L dcj) L N B 



N 2 T ds r 



N ds r , 



N dSr, 



(B12) 



Thus, by consideration of the elliptic integral representa- 
tion, it can be shown that 



1 dT 

JPTdVr, 



1-sj 



(B14) 



where 



^l/r _ dcf> L / R ds g d(/) L/R 



ds r , 



ds r 



+ 



ds„ ds ri 



and from this result follows the force on the pin given in 
(B13) Eq. pig)). 
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